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We study electron transport through multiply connected mesoscopic geometries containing interact- 
ing quantum dots. Our formulation covers both equilibrium and nonequilibrium physics. We discuss 
the relation of coherent transport channels through the quantum dot to flux-sensitive Aharonov- 
Bohm oscillations in the total conductance of the device. Contributions to transport in first and 
second order in the intrinsic linewidth of the dot levels are addressed in detail. We predict an 
interaction-induced asymmetry in the amplitude of the interference signal around resonance peaks 
as a consequence of incoherence associated with spin-flip processes. This asymmetry can be used to 
probe the total spin of the quantum dot. Such a probe requires less stringent experimental conditions 
than the Kondo effect, which provides the same information. We show that first-order contributions 
can be partially or even fully coherent. This contrasts with the sequential-tunneling picture, which 
describes first-order transport as a sequence of incoherent tunneling processes. We predict bias- 
voltage-induced Aharonov-Bohm oscillations of physical quantities which are independent of flux in 
the linear-response regime. Going beyond the Onsager relations we analyze the relations between 
the space symmetry group of the setup and the flux-dependent nonlinear conductance. 



PACS numbers: 73.23.Hk, 73.63.Kv, 73.40. Gk 
I. INTRODUCTION 

The manifestation of quantum coherence in finite sys- 
tems is in the foundations of the physics of mesoscopic 
systems. The presence of quantum coherence is de- 
tectable through interference experiments, most notably 
Aharonov-Bohm (AB) interferometry. The prototype of 
an AB setup is a double-slit experiment, as shown in 
Fig. |lja. An electron moving from the left reservoir to 
the right one is split into two partial waves which inter- 
fere with each other. A magnetic flux $ that penetrates 
the area enclosed by the two paths changes the relative 
phase of the amplitudes of the two partial waves, which 
yields a flux dependence of the total transmission prob- 
ability through the device. 

A convenient framework to study the role of electron- 
electron (e-e) interactions in mesoscopic systems is sug- 
gested by employing quantum dots (QDs). The lat- 
ter provide a relatively simple and controlled scheme to 
address e-e interactions without undermining the rich 
physics involved. Transport through QDs has been stud- 
ied extensively and revealed interesting phenomena such 
as resonant tunneling, Coulomb blockade, and the Kondo 
effect. However, the measurement of the current through 
QDs does not provide information about the quantum 
coherence of the transport. In particular, it doesn't ad- 
dress the interplay of electron-electron interactions and 
quantum coherence. In order to approach these qWQSn 
tions QDs have been embedded in AB geometries.ErEHl 
Magnetic-flux sensitivity of the total current has been 



observed ,@flOL3'E3 indicating at least partially coherent 
transport through the QD. 

By applying a finite bias voltage across the device the 
QD can be driven out of equilibrium, lifting the complex- 
ity of the system to a qualitatively higher level. In the 
context of interacting QDs nonequilibrium effects such 
as splitting of Kondo resonances have been discussed. It 
is only natural to expect that driving the system out of 
equilibrium enriches the interplay between quantum co- 
herence and interaction effects by another facet. Some of 
the symmetries present at equilibrium, which underline 
linear response, may be broken, and at the same time 
new qualitative features may emerge. 

In this paper we address two main issues: 

(i) What can we learn about QDs by embedding them 
into AB geometries? Based on Ref. Ill], we will classify 
different contributions to linear transport into coherent 
and incoherent channels. We concentrate on the domi- 
nant transport mechanism, which is of first order in the 
intrinsic linewidth of the dot levels when the QD is tuned 
in resonance with the leads and of second peder ( "cotun- 
neling") in the Coulomb blockade regimeO A sufficient 
condition to establish full coherence is to find an AB ge- 
ometry in which the interference pattern reveals fully de- 
structive interference since transport processes through 
the QD which are cancelled by adding parallel transmis- 
sion channels cannot be incoherent. As we will see the 
proper choice of the AB interferometer, including cither 
one or two QDs, is important for the manifestation of 
this criterion. 

A very important and striking result of our analysis 



will lead to the suggestion that AB interferometers (with 
either one or two QDs) can be used to probe the total spin 
of the QD. We will find that the relative magnitude of 
the flux-dcpcndcnt AB signal, the visibility, changes from 
one Coulomb blockade valley to the next. This amounts 
to an asymmetry in the AB signal around the Coulomb 
conductance peaks. The higher visibility indicates the 
side of the peak which corresponds to a situation where 
the QD is mostly occupied by an even number of elec- 
trons (doubly occupied levels) with a total spin 0. The 
lower visibility corresponds to an odd number of elec- 
trons and total spin 1/2. A similar even/odd effect in the 
low-tempcrature,e«]jiductance through a QD results from 
the Kondo effect B3iH3 Thus, measurement of the conduc- 
tance probes the total spin of the QD, too, although the 
physical origin is completely different. We stress that our 
spin-dependent asymmetry effect is not a manifestation 
of incipient Kondo physics. Indeed, we expect our asym- 
metry to disappear in the zero-temperature limit, where 
the Kondo effect is fully developed. In order to access 
the Kondo regime, low temperature and strong coupling 
of the dot to the leads is required. This contrasts with 
our suggestion to employ an AB interferometer. We pre- 
dict the asymmetry of the AB signal to survive at high 
temperature and weak coupling, i.e., much less stringent 
conditions. 

(ii) The other major issue of this paper addresses 
the interplay between quantum coherence and electron- 
electron interaction for systems out of equilibrium. In 
the linear-response regime the average electron number 
on the QD is an equilibrium quantity. For weak coupling 
between dot and leads (such that the intrinsic linewidth 
is small compared to the energy scale provided by the 
temperature and the level spacing), the dot occupation 
is governed by classical Boltzmann factors, which are in- 
dependent of magnetic flux. Out of equilibrium, however, 
we predict bias-voltage-induced AB oscillations in the dot 
occupation for a single-dot AB interferometer. We em- 
phasize that, again, the experimental conditions for this 
effect to occur, high temperature and weak coupling, are 
accessible rather easily. 

Qualitative differences between equilibrium and 
nonequilibrium situations also show up in symmetry re- 
lations of transport coefficients. The classical double-slit 
setup (Fig. nk) is an example of an open- geometry inter- 
ferometer, where those electrons which are not transmit- 
ted directly to the lead on the right are absorbed by other 
gates and terminals of the system at the periphery of the 
device (these are not shown in the figure). In most parts 
of the paper, however, we consider AB interferometers 
in a closed geometry, as sketched in Fig. nb, where these 
electrons can only go to either the terminal on the left or 
the one on the right-hand-side. For the open-geometry 
setup (Fig. [j]a) <f> = is not a symmetry point (i.e., the 
transmission, as well as transport coefficients and the 
partition function, are not invariant under <f> — > — $). 
The symmetry point depends on specifics of the interfer- 
ometer. This contrasts with the case of a two-terminal 



geometry (Fig. |Tjt>) , where all thermodynamic potentials 
as well as linear-response transport coefficients are in- 
variant under the inversion of the flux.c3 The fact that 
$ = is an extremum point is usually referred to as phase 
locking, and is a direct consequence of Onsager relations. 
It has been demonstratecu that beyond linear response, 
i.e., in the presence of a finite voltage bias, phase locking 
is, in general, broken. We show by an explicit calcula- 
tion that for a two-terminal AB interferometer including 
a single QD both electron-electron interaction and finite 
bias voltage is needed to break phase locking. 

But even in the nonlinear-response regime, there are 
symmetry relations in the transport coefficients which 
are connected to the space-symmetry group of the setup. 
This will be especially important for AB interferometers 
with two QDs, where, depending on the spatial symme- 
try, phase locking can occur even at finite bias voltages. 
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FIG. 1. a) Two-slit experiment as an open-geometry AB in- 
terferometer: only a fraction of the electrons passing through 
the interferometer will reach the drain, b) Two-terminal 
(closed-geometry) AB interferometer: all electrons injected 
from the source must either reach the drain or be reflected to 
the source. 

We conclude this introduction with some remarks 
about the formalism to be employed. The experimen- 
tally accessible quantity which characterizes transport is 
the total current through the device, which is driven by 
bias voltage V. The electrical current (flowing from the 
right to the left) can be written in the form 



'-££ 



dwT <r («)[/ £ (w)-/,t(a;)], (1.1) 



where / r (w) = l/[l+exp(/3(w — H r ))] with r = L/R is the 
Fermi-Dirac distribution function, and (3 — l/(fcsT) is 
the inverse temperature. The left and the right lead have 



fixed electrochemical potentials (j,l = /-i + eV/2 and [/,r = 
(i — eV/2 for the electrons, respectively, with e > 0. The 
quantity T a {u>) defines the probability for an incoming 
electron with energy u) (measured with respect to fj,) and 
spin a to be transmitted through the device. 

Accordi ng , t o the Landauer-Biittiker or scattering 
formalism,? 3 HP 2 ! which we are not going to use for rea- 
sons explained below, the transmission probability T a {ui) 
for each channel labeled by u> and a is obtained from the 
modulus squared of the total transmission amplitude. In 
the example of the double-slit setup, Fig. Ilia,, we have 
T a (uj) — \tia(u)) + t 2c r{uj)\ 2 , where t\ and t 2 are the 
partial amplitudes through cither slit and are complex 
quantities. The associated phases have an orbital contri- 
bution, determined by the geometrical details of the in- 
terferometer, and a magnetic-flux-dependent part. The 
transmission probability is, then, given by 



T=\h\ 2 + \t 2 \ 2 +2\t 1 h\coa(<p + 50) 



(1.2) 



where 69 is the relative orbital phase, ip = 27r$/$o ac- 
counts for the enclosed magnet flux, and <I>o = h/e is 
the flux quantum. (Here, for this geometry, it is jus- 
tified to neglect the multiply scattered higher-winding- 
numbcr amplitudes. For closed geometries, such as in 
Fig. [l]b, the partial amplitudes may exhibit a more com- 
plicated dependence on magnetic field, representing par- 
tial amplitudes with higher winding numbers around the 
enclosed flux.) The success of the Landauer-Buttiker pic- 
ture lies in part in the fact that, once the partial trans- 
mission amplitudes are known, the analysis is straightfor- 
ward. Indeed, for noninteracting systems, the transmis- 
sion amplitude for an electron to travel from the left to 
the right lead is simply determined from the (retarded) 
single-particle Green's function G r LRa {u>) which involves 
a creation operator of an electron in the right and an an- 
nihilation operator in the left lead. This single-particle 
picture, however, breaks down if electron- electron inter- 
action is taken into account. This has been shown in the 
literatures and will be explicitly demonstrated in the 
next section. In this case, we do not see a general recipe 
of how to construct meaningful transmission amplitudes. 

Instead we directly calculate the current usiagJ^xeen's 
function techniques for interacting systems J^'lj'Lj and 
extract the transmission probability from Eq. fll.ip with- 
out worrying about transmission amplitudes. 

The outline of this paper is as follows. In Sec. O we 
employ a simple physical picture to distinguish coherent 
from incoherent transport through a QD. To substantiate 



III 



these qualitative considerations we calculate in Sec 
the current through an AB interferometer in which a sin- 
gle QD in embedded. Interferometry with two QDs is 
addressed in Sec. [V. Sectio n yj\ is devoted to nonequi- 
librium physics, and in Sec. VI we draw the connection 



between the spatial symmetry of the setup and symme- 
try relations of nonlinear transport coefficients. We have 
also included a few appendices to provide some technical 
details of our analysis. In Appendix |A| and [B| we deter- 
mine the density matrix for the single-dot and two-dot 



AB interferometer, respectively, and in Appendix y we 
present an alternative derivation of Eq. (4.7). 

The m ain result s of o ur analysis are summarized in 
(i) Eqs. (3.9) and ( 3.12 ) (the zeroth and first-order ex- 
pressions for t he cur rent through the interferometer); 
(ii) Eqs. fl3.13| ) and ( 3.14 ) (the flux-dependent linear- 
response transmission for the clo sed a nd op en geo metry 
respectively); (iii) Eqs. fljUJ ), J3T9J) , Q , Q (the 
absence/presence of asymmetry in the case of U = 0, 
U = oo respectively), for the single QD in the cotun- 
neling regime, and the double QD at resonance; (iv) the 
various regimes for the multilevel dot indicated at the end 
of Section |IIIE| ; (v) Eqs. (|1|) and (§!]) (the finite-bias 
induc ed A B osc illations in the dot occupation (N)); (vi) 
Eqs. (6.5) - (3/7) (the relation between flux-reversal sym- 
metry and space-symmetry, away from linear response) . 



II. COHERENT VERSUS INCOHERENT 
TRANSPORT THROUGH A QD 

Before analyzing AB interferometers which contain 
QDs, as shown in Figs. 0b and c, we begin with con- 
sidering transport through a QD alone, Fig. ga. 
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FIG. 2. Transport through a QD and AB interferometers 
with either one or two QDs. 



We develop a simple physical picture to underline the 
essential physics involved in the predicted interaction- 
induced asymmetry of the interference signal around res- 
onance peaks. All statements made in this section will be 
backed up by strict calculations for the AB interferom- 
eters in the subsequent sections. Furthermore, we give 
an explicit example of how the single-particle formalism 
breaks down in the presence of Coulomb interaction. 

We consider a single-level QD with level energy e, mea- 
sured from the Fermi energy of the leads. The Hamilto- 
nian 



rdot 



H am =H l +Hr + H d + Ht 



(2.1) 



contains a part describing noninteracting electrons in the 
left and right leads, 



U r — 



E 

ko 



t-kr dfco-r ^kcrr 



with r = L/R. The isolated dot is described by 



H D = eJ2', 



Un-fni 



(2.2) 



(2.3) 



where n a — c\c a counts the number of electrons with 
spin a. The energy of the dot level e can be varied by an 
applied gate voltage. The electron-electron interaction 
is accounted for by the charging energy penalty U — 
2(e 2 /(2C)) for double occupancy, where C is the effective 
capacitance of the QD. To keep the discussion simple we 
choose the generic limits U — for the noninteracting 
case and U = oo for an interacting QD. We stress that 
the pursuant analysis, and in particular our qualitative 
conclusions, are applicable to the case of finite U as well. 
Tunneling between the QD and the leads is modeled by 



Ht = y> r o£„.C <r + H.C. 



kcrr 



(2.4) 



(We neglect the energy dependence of the tunnel matrix 
elements t^/RJ. Due to tunneling the dot level acquires 
a finite linewidth T = T^ + Y r with 



L/R 



MtL/ R \ 2 N L/R 



(2.5) 



where N L i R is the density of states in the corresponding 
lead. 

For the sake of developing physical insight, we consider 
an off-resonance scenario, i.e., the gate voltage is tuned 
in such a way that |e| 3> fcgT, T. In this regime (and 
as long as the temperature is higher than the Kondo 
temperature), transport is dominated by second- 
processes in r, which are usually called cotunnelingi 
Furthermore, we assume an infinite large charging en- 
ergy U — oo. Figure || shows all processes (second order 
in T) which contribute to the current from left to right 
in which the incoming electron has spin up. Similarly, 
there are processes in which the incoming electron has 



spin down. They will introduce a trivial factor of 2 to 
the final expressions for the current. The three possible 
second-order processes are 

(a) a (spin-up) electron enters the QD, leading to a vir- 
tual occupancy, and then leaves it to the other side. 

(b) a (spin-up) electron leaves the QD, and an electron 
with the same spin enters. 

(c) a (spin-down) electron leaves the QD, and an electron 
with opposite spin enters. 



a) 





b) 



c) 




FIG. 3. Cotunneling processes for U = oo. The solid line 
indicates the process that happens first, the dashed line the 
process that occurs afterwards. Double occupancy of the QD 
in the initial, intermediate, or final state is prohibited by 
charging energy. 

Note that double occupancy of the dot (even in a vir- 
tual state) is forbidden since we have assumed U = oo. 
All processes are elastic in the sense that the energy of the 
QD has not changed between its initial and final states. 
Process (c), though, is incoherent since the spin in the 
the QD has been flipped. In other words, the travers- 
ing electron has left a trace in the "environment" (i.e., 
the dot). This implies that when the QD is embedded 
in one of the arms of an AB interferometer, one will be 
able to tell that an electron participating in process (c) 
has indeed traversed the arm with the QD (rather than 
the other arm). Such a process will then contribute to 
the total current but not to the flux-sensitive component 
thereof, independent of the realization of the AB inter- 
ferometer. The notion, that energy exchange is not nec- 
essary for dcphasing,E3 and that the latter can take place 
through, e.g., a spin flip of as, external degree of free- 
dom, has been made early on.Ea In our case the electrons 
in the |Qp itself (and their spin) serve as the "dephasing 
bath" .El 

Our simple picture already indicates that the physics 



in adjacent Coulomb blockade valleys is qualitatively dif- 
ferent. This is the origin of the asymmetry in the visibil- 
ity discussed throughout this paper. In our model there 
are two possible charge states for the QD: the dot being 
empty, N = 0, or singly occupied, N = 1, with cither spin 
up or down. By tuning the gate voltage (which changes 
the level position e) one can drive the system away from 
resonance, where both states N — and N = 1 are possi- 
ble, to the Coulomb blockade valley with a fixed electron 
number. The latter is either N — or, on the other 
side of the resonance, N — 1. In the first case, the main 
contribution to the current is due to the coherent pro- 
cess (a) in Fig. ||. In the latter case, two processes, (b) 
and (c), contribute to the current but only one of them, 
namely (b), is coherent. We thus expect that in an AB 
interferometer with a single QD the AB amplitude at 

a gate voltage V g in the N = valley will be twice 
as large as the AB amplitude at a gate voltage Vg in 
the N — 1 valley. For such a comparison one needs to 



,(0) 



(i) 



consider gate voltages V g and V g for which the flux- 
averaged transmission is the same. We also note that a 
multilevel system may offer more complex configurations, 
e.g., a resonance which separates a valley with all levels 
doubly occupied from a valley where one extra electron 
is added and one, three, or more levels are now singly 
occupied. 

The transmission probabilities of electrons with energy 
uj near the Fermi level of the leads can be obtained by 
calculating the transition rate in second-order pertur- 
bation theory and multiplying it with the probabilities 
P x to find the system in the corresponding initial state 
X- For an incoming electron with spin up the trans- 
mission probabilities are P X T lT nRe[l / (ui — e + i0 + ) 2 ] 
with x = 0, t,4 for case (a), (b), and (c), respectively. 
Since P + P T + P[ = 1 and P + P a = 1/[1 + /(e)] in 
equilibrium, where /(e) is the Fermi function, we find 
for the transmission probability throughJJie QD that 
l£ ot (w) = T dot > coh (w) + T^ ot ' incoh (u;) withEl 



T dot M 

Tidot,coh/, a 



Re 



r r ,r 



Li- R 



(uj — e 

j; dot n 

i + /W 



M+y 



(2.6) 
(2.7) 



In linear response the participating electrons have ener- 
gies uj that are spread by fc^T around the Fermi energy. 
Since we consid er a n off-r eson ance situation, we can set 
U) ~ in Eqs. (21) and (|j) and find that, while the 
total transmission is symmetric under e — > — e, the factor 
1/[1 + /(e)] ascribes an asymmetry to the coherent part. 
We have argued that incoherence is induced by flipping 
the spin of the transfered electron, Fig. ||c. One might ex- 
pect naively that such spin-flip processes may take place 
for noninteracting systems as well. On the other hand, 
in the absence of interaction, transport should be fully 
coherent. This puzzle is solved by the observation that 
in the absence of interaction, U = 0, double occupancy 



of the dot is allowed and, therefore, more processes are 
possible. Figure [| depicts these additional second-order 
processes (with an incoming spin-up electron). Spin flip 
takes place in Figs. |^c and [Ip only. The first is hole-like 
while the second is particle-like (the solid line indicates 
the first step and the dashed line the second step of the 
process) . Both amplitudes have the same magnitude but 
come with opposite signs and, hence, cancel out. 



a) 



b) 



c) 




FIG. 4. Additional cotunneling processes for U < oo. Now 
also double occupancy of the QD is possible. 

Thus, in the absence of Coulomb interaction (and no 
coupling to an external environment) there are no spin- 
flip processes, and the transmission is fully coherent. We 
perform an analogous second-order perturbation calcula- 
tion as above, use Po + Pp + Pj_ + P-p^ = 1. and find that 
T^ ot (w) = T^ ot ' coh (ui) given by Eq. (|1|) without any 
asymmetry factor. Once U is turned on, the exact can- 
cellation of the spin-flip processes is lost, and incoherent 
transmission channels induce an asymmetry of the AB 
signal. 

As mentioned in the previous section, for noninter- 
acting systems the transmission amplitude is determined 
by the retarded single-particle Green's function from the 
left to the right lead, and then we can easily apply the 
Landauer-Biittiker formalism. In our case this amplitude 
can be expressed as 



tf t {u) = i^T L Y R G r a {u), 



(2.8) 



where the dot Green's function G r a (ui) is the Fourier 
transform of — i<d(t)({c a (t), c^(0)}). We insert the dot 
Green's function (for U — and zeroth order in T) into 
Eq. (|J) and get 4 ot (uj) = i/LYF^/(w - e + i0+) which 
yields 



dot/ 



ie>)r=7r», 



(2.9) 



with T^ ot (w) given by Eq. (gj). 

In £/ie presence of interaction, however, Eq. (2.L) is no 



longer a good definition for the transmission amplitude. 
To see this we calculate the dot Green's function in the 
limit U — oo and zeroth order in T. We find t^ ot (w) = 
i(P + P T )v / Tlr^/(w - e + i0 + ) and, therefore, 



ICW = 



[l + /(e)] s 



(2.10) 



We see that not onl y do es |i£ (w)| not yield the to 



tal transmission Eq. (2.6) through the dot, but i t di ffers 
from the coherent part of the transmission, Eq. (2.7), as 
well: there is no direct physical meaning of the expression 

Can we find a similar clear picture for first-order trans- 
port in r, which dominates at resonance? As long as 
quantum coherence is not addressed, first-order trans- 
port through QDs has been successfully described by the 
use of the sequential-tunneling picture or "orthodox the- 
ory" . Tunneling rates into and out of the QD are cal- 
culated with Fermi's golden rule. They are connected to 
the probabilities for the dot states by a master equation, 
and determine the current. According to this philosophy, 
different tunneling processes in or out of the QD are un- 
correlated since it is implicitly assumed that the electron 
experiences some dephasing during its stay in the QD. 
As a consequence the sequential-tunneling picture implies 
that first-order transport should be completely incoher- 
ent. This may be well justified for large QDs which are 
coupled to phonons or other degrees of freedom. But for 
the system under consideration in this paper, a single- 
level QD with no coupling to external baths, this is not 
the case. For U = we even expect fully coherent trans- 

As we will prove in Sec. IV by 



port to all orders in T. 
strict calculations, the sequential-tunneling picture com- 
pletely misjudges the coherence of first-order transport, 
although it predicts the correct transmission probability 
through a QD. In our case, first-order transport should 
rather be viewed as "resonant tunneling" expanded up to 
first order than visualized as a sequence of (incoherent) 
sequential-tunneling processes. 



III. INTERFEROMETRY WITH A SINGLE QD 

To support the results of our intuitive picture, we an- 
alyze quantitatively AB interferometers which contain a 
single QD, such as shown in Fig. pp. We emphasize that 
all following resul ts (except for the discussion of many 
levels in Sec. Ill Ej ) are derived by strict calculation of the 
total current through the device in the presence of mag- 
netic flux. At no point we make use of the qualitative 
picture outlined in the previous section. The compati- 
bility of our calculations with the intuitive picture serves 
rather as a check. 

Electrons emitted from the left lead have two possible 
ways to reach the drain on the right. They can either 



go through the QD in the lower arm or choose the upper 
arm. We denote the transmission amplitude for the latter 
trajectory by i rof and assume that it is independent of 
energy u> and spin a. The total transmission probability 
T^ ot (cj) through this device is the sum of three parts, 



Tl ot {uo) = T« ot {uj) + T™ + T* ux (uj) . 



(3.1) 



The terms T^ ot (w) and T* ci = \t lci \ 2 are the flux- 
insensitive transmissions through the QD and the ref- 
erence arm, respectively. (In principle, T^ ot (uj) feels the 
presence of the reference arm and vice versa. When low- 
order transmission in the couplings is considered this in- 
fluence may be ignored.) Interference is described by the 
remaining term, T^ ux (oj), which depends on magnetic 
flux. 



A. Current formula 

In the first step we relate the flux-dependent nonlinear 
current to the dot Green's function for the closed geom- 
etry, Fig. gb. This derivation is somewhat technical but 
straightforward. The central res ult of this part is the 
current in first order in i rof , 
zeroth-order result, Eqs. (3.9)]. 



Eq. ( |3.12| ) [along with the 
The total Hamiltonian, 



H = H 



dot 



H 



ref 



(3.2) 



consists of two parts. The first, H dot , given by Eq. ( pj 
describes the arm containing the QD, and the second, 



H 1 



ref 



E (*° 



kaR L 



qaL 



H.c.) 



(3.3) 



k£R,qEL,<7 



with 2iri^JNLNR — It'^e 29 , models the transmission 
through the reference arm. In general, the magnetic flux 
$ enters the phases of all three tunneling matrix ele- 
ments £l, t R , and t Tcl . Above, however, we have chosen 
a gauge in which only t acquires a flux-dependent phase 
ip = 27r$/$o but leaves tr, and t R flux independent, and 
we can choose the latter to be real. 

The operator for the current from the right lead is 
given by the time derivative of the total electron number 
operator n R = ^2 keR a a \ aR a kaR times the elementary 
charge e. This yields for the total current 



I = Ir = e = i-([H,n R \ 

at n 



(3.4) 



The latter expression yields Green's functions which in- 
volve Fermi operators of the right lead and of either the 
left lead or the dot, 



'« = -! E 



q£L,keR,rr 



djj 



iG<(cu)+R.c. 



keR,a 



~\ E dcj\t R G< k Ju;)+U. C . 



(3.5) 



with the notations Gf k a (t) 



,t 



*«rn(0)<Vi(*)) and 7 (o) = 



The indices g and fc label 
the states in the left and right leads, respectively. The 
index d indicates that a dot electron operator is involved 
(in our simple model there is only one dot level). The 
other Green's functions are defined similarly, G k a (t) = 

*(oU(°)°*'*(*)) and G kdJ f ) = *(4(0)o tofl (*)); Note 
that, since the z component of the spin of the entire sys- 
tem is conserved (as well as the total spin), only Green's 
functions which are diagonal in spin space are involved. 

The first (second) line of Eq. (pit) describes electron 
transfer from the left (from the QD) to the right lead 
or vice versa. This transfer can be a direct tunneling 
process or a complex trajector y th rough the entire de- 
vice. We emphasize that Eq. ( |3.5| ) is exact as long as 
the full Green's functions with contributions of arbitrary 
high order in t and t are inserted. 

Our goal is to derive a relation between the current 
and Green's functions involving only dot operators. To 
achieve this we employ the Keldysh technique and use 
the matrix representation 



G r G< 
G a 



(3.6) 



where G r and G° are the usual retarded and advanced 
Green's functions, respectively. 

The contribution to the total current can be classified 
by powers of t Tct . The zeroth-order term I R ' represents 
the current through the QD in the absence of the refer- 
ence arm and is independent of magnetic flux. The main 
contribution I R to the flux-dependent part is first order 
in t . In the spirit of a series expansion we drop in the 
flux-dependent part all second- or higher-order terms of 
t rcl , which is a good approximation as long as the trans- 
mission through the reference arm is small (trajectories 
with higher winding numbers, where an electron moves 
around the enclosed flux several times, are described by 
such higher-order terms). Transmission through the ref- 
erence arm, which is at least of order (i rcf ) 2 , is not con- 
sidered in this subsection. 



We start with the zeroth-order term, I 



(0) 

R 



Only 

the second line of Eq. (3.5) has to be included since 
the first line explicitly contains i. Since the electrons 
in the leads are noninteracting, we get the Dyson-like 
equation G^ >0 .(w) = Gh °\u))t R gk,a{u>), where G<(t) = 
i{c\c a {t)) is the Green's functions of the dot, and g ka = 

*( a fcCT_R a fe<rfl(^)) f° r tnc leads, and the retarded and ad- 
vanced Green's functions are defined similarly. We do 

the analog for G fed . For the (noninteracting) leads we 

make use of g kiCT (uj) = 2nif R (uj)8(uj - e k ), g^ a (oj) = 

l/(w - e k + i0+), and g-Juj) = (<£> 



>£ 



du> 



G <(o) 



h. (g: 



(0) _ Qa(0) 



(3.7) 



As a check we verify that in equilibrium, V = 0, the 
relation G < (u) + f(uj)[G r (uj) - G a (uj)} = guarantees 
that no current is flowing, I R = 0. 

It is straightforward to derive an analogous expression 
for I r for the left lead. The total current must be con- 



served, I L 



r(0) 

L R 



0. Even more, the current must 



be conserved for each spin separately. This yields the 
condition 



0= / dtu 



Q<(0) 



Tl/l + T R f R 
Tr.+Tn 



/-ir(0) _ fm(0) 



(3.8) 



for the integral (but not for each energy uj separately). 
Due to current conservation we can write I R = (T lI R ~ 
T r Il)/(Tl + Tr) and find the well-known result for the 
current through a QD (in the absence of a reference arm) , 



r(0) _ _ 2e 

1 R — 



T,V 



Ll_R 



h Yi 



YsJdojImG^ih-fR) 



(3.9) 



and A a (uj) = — (1/ir) ImG^(w) is the spectral density. 

We now turn to the flux-dependent term I R . Similarly 
to the above, we can make use of Dyson-like equati ons . 
We get a contribution from the first line of Eq. (3.5) 
and use gJJ„ v 



i(0), 



r (w) = gq,cr(v)tLGlr (u>)t R gk,v(u>). For 
the contributions of the second line of Eq. ( |3.5| ) we 
make use of G^ (T (o;) = Gi 0) (u))t L g qt<J (u))ig k ,cr(u) + 

Ga (u>)tRgk,cr(w). After collecting all terms and em- 
ploying-.the relation ±it 2 it R t L N L N R = ^T L T R \t loi \e iip 
we getc3 

& = Zy/r^\r t \ , Z J du, 

'coa^(G^+G^)(f L -f R ) 
+zsin^(G<( )+/ L (G;(°)-G^)))" 

-f r «E/ d - [G< W +f R (G^-G^)} . (3.10) 

Again, we can check that due to the equilibrium re- 
lation G<(u) + f{uj)[G r {uj) - G a (ui)} = the current 
vanishes at zero bias voltage. The current conservation 
(for each spin) in first order in t lcl is equivalent to 



vFlFr , f| 



We obtain then 



r L + r R 

du) 



\f ct \ sin^ J dLo (CTJ® - Gf>) (f L - Jr) 



r <(i) T L f L + TrJr 



(3.11) 



Note that to ensure current conservation, Green's func- 
tions of higher order in t rcf are important since the first 



line of Eq. (|3.1l|) alone is in general nonzero (compari- 
son with Eq. fl3.9| ) shows that it is even proportional to 
the current I^ in zeroth order). This means that for 
a systematic and consistent description the influence of 
the reference arm on the QD has to be taken into ac- 
count properly. It is, in general, not sufficient to treat 
the Green's function of the QD independent of the sur- 
rounding environment in the device. One has, therefore, 
to check carefully in each case up to what extent the 
interference signal probes the properties of the QD. 

We use I R = (r L I R -T R I L )/{r L +T R ) and the current 
conservation in zeroth order, Eq. ( |3.8| ), to get 

Ir = ^VrZf^|i ref | cos^ Wdu, ReG;(°)(/ L - f R ) 

-|r^ff7^/^ Im ^ (1)(/L - /fl) - (3 - 12) 

a 

This is the central result of the present derivation. Again, 
we emphasize that the influence of the reference arm on 
the QD plays a role (it shows up in the Green's functions 
which are of first order in t ref ). Furthermore, we remark 
that the Green's functions G^ and G^ include multiple 
tunneling between QD and leads, i.e., they include con- 
tributions from arbitrary high order in t (the superscripts 
(0) and (1) only labels the order in t ). 

B. Linear- response regime 



In the linear-response re gime , we replace /l(^) — 
f R {uj) -* -eVf'{u) in Eq. (|3.12f ) and take the G reen' s 
functions at equilibrium. The first line in Eq. (3.12) 
is obviously symmetric under reversal of magnetic flux, 
ip — » —ip. Moreover, in equilibrium also the Green's func- 
tions G^ have this property.e3 This establishes phase 
locking in linear response, as expected for a two-terminal 
device. 

The situation simplifies further if we concentrate on 
the dominant contribution to the flux-dependent linear 
conductance, which is of first order in T and first order 
of . In this case the second line of Eq. (3.12) drops out 



in f 



(since there is a prefactor of order T and the equilibrium 
Green's function G' 1 ' in zeroth order in T vanishes). Th e 
connection to the transmission is established by Eq. (1.1). 
We obtain 



T CT flux H = 2^/T L T R ~\f ct \ cos^ReG^M , (3.13) 

and realize that up to the given order the influence of 
the reference arm on the QD does not matter for the 
transport properties. 

For comparison we cite the result for a system with an 
open geometryJU 



with 9 = ip + A9, where A6 is determined by the sp ecific s 
of th e int erferometer. The two expressions Eqs. ( 3.13| ) 
and (3.14) look very similar. One should keep in mind, 



though, that in higher orders in T or in nonlinear response 
an add itiona l term, which originates from the secon d line 
of Eq. ( 3.12J ) enters, while for the derivation of Eq. ( 3.14 ) 
one considers a special case of open geometry where the 
reference arm has no effect on the QD's properties. 

The transmission Eq. ( |3 .13 ) is always extremal at 
if = 0. Such a "phase locking" does not take place in 
the open-geometry setup: the AB phase at which the 
transmission is extremal can be continuously varied by 
tuning the energy of the dot level via a gate electrode. 



C. Probing the total spin of an interacting QD 

In linear response and lowest order in t and F we 

get G r J 0) (Lu) = l/(w - e + i0+) for U = and, accord- 



ing to Eq. (3.13), the amplitude of the flux-dependent 
linear conductance is symmetric around e = 0. This is 
consistent with our intuitive picture, since in the absence 
of spin-flip processes no incoherent channels are involved 
on either side of the resonance. 

This changes when electron-electron interaction in the 
QD is important. For U — oo the retarded Green's func- 
tion 



g; (0) ( 



i 



u 



l + /(e) uj-e + i0+ 



(3.15) 



acquires a prefactor P + P a = 1/[1 + /(e)] in linear 
response and lowest order in i rcf and T. This prefactor 
indicates an "interaction-induced" asymmetry associated 
with spin-flip processes. 

For illustration we show in Fig. H the flux-dependent 
linear conductance as a function of the level position e 
for a fixed value ip = 0. In order to obtain the flux de- 
pendence, the result simply has to be multiplied by cos ip. 
At e = the AB oscillations vanish. [ Note that we con- 
sider here only the first harmonic Eq. ( |3.12 ). Higher har- 
monics survive, an effect known as frequency doubling.] 
In the absence of interaction, [7 = 0, the amplitude of 
the AB signal is symmetric around that point, while for 
U — oo clear asymmetry is predicted. When the QD is 
most likely filled up with one electron, the total spin of 
the QD is 1/2 and spin-flip processes reduce the interfer- 
ence signal. On the other side of the resonance, e > 0, 
the QD is most likely empty. The total spin is and the 
AB oscillations are by a factor of 2 larger than for e < 0. 



T, 



flux, open 



H = 2VTlT^|i ref |Re e- iB GrW(uj) , (3.14) 



u = o 

infinte U 




dot 



FIG. 5. Interference signal for AB interferometer with a 
single QD. Asymmetry in the magnitude of the signal appears 
for nonzero U. Plotted is the flux-dependent part of the con- 
ductance, normalized by i ref |r/(fcsT), in units of e 2 /h. 

It is well known that due to the Kondo effect the trans- 
port through a QD increases at low temperature for e < 
but not for e > 0. Therefore, the Kondo effect reveals the 
same information as our procedure. We note, however, 
that the Kondo effect occurs under much more stringent 
experimental conditions, namely, strong coupling of the 
QD to the leads and low temperatures (in comparison 
to the Kondo temperature). In contrast, our approach 
to detect the spin of the QD works in the experimen- 
tally easier accessible regime of weak coupling and high 
temperatures. 



D. Coherence of non-spin-fiip cotunneling 

As has been outlined in the introduction we use the fol- 
lowing procedure to investigate the coherence of cotun- 
neling. We tune the system away from resonance, where 
transport though the QD is dominated by the usual co- 
tunneling processes. We then try to find an AB setup 
such that the total conductance vanishes once a refer- 
ence arm is added, i.e., we look for complete destructive 
interference of the coherent channels. (The existence of 
AB oscillations alone only proves partial coherence.) We 
remind ourselves that the transport through the AB in- 
terferometer probes many channels, characterized by en- 
ergy to and spin er, simultaneously. Zero transport can 
only be achieved when all channels show destructive in- 
terference at the same time. To achieve this we want to 
adjust the magnitude of the amplitude t for transmis- 
sion through the reference arm such that T^ ot (uj) = T^ ci 
for all contributing energies. 

Cotunneling is the dominant transport channel when 
|e| ^$> r, ksT applies. In this regime, the linear-response 
conductance through the QD in the absence of the refer- 
ence arm is 



0I__ 
dV 



2e 2 T L T R 



(3.16) 



v=o 



for both U = and U — oo. Cotunneling through the 
QD is of second order in T. The conductance through 
the reference arm in the absence of the QD is 



dr ci 



dV 



2e 2 



,ref|2 



(3.17) 



v=o 



We now make the adjustment |£ | = ^/Yi,Tr/\(\, i.e., 
\t lcl | is of order T. Therefore, the Green's function which 



enters the flux-dependent transmission Eq. (3.13) is of 



zeroth order in T and t , i.e., the influence of the ref- 
erence arm on the QD is not probed under the present 
conditions. We find for the total conductance 



Qrtot 



dV 



4e 2 r L r, 



v=o 



1 - -r-r cos <p 
e 



for the noninteracting case and 



dP' 



dV 



= 4- 



T,Y 



L 1 R 



V=0 



COS^ 



=| 1 + /(e) 



(3.18) 



(3.19) 



for U — oo. This shows that cotunneling in the noninter- 
acting case is fully coherent. In the interacting case spin- 
flip processes are present which spoils coherence. This is 
described by the asymmetry factor l/[l + /(e)], in accor- 
dance with our intuitive picture. 

We conclude that an AB interferometer containing a 
single QD is suitable to prove coherence of non-spin-flip 
cotunneling through the QD. 



E. Many levels 

In the model we discussed so far, the spacing A of the 
dot levels was assumed to be larger than the charging en- 
ergy and the energy scale provided by the temperature, 
the intrinsic linewidth or the voltage bias, so that only 
one dot level participates in transport. We found that 
the ratio of the number of coherent to the total num- 
ber of cotunneling channels was 1 or 1/2 in the valley 
where the electron number of the dot is even or odd, re- 
spectively, giving rise to an observable asymmetry in the 
interference signal. The resulting sequence of asymmet- 
ric AB oscillations is shown schematically in Fig. g. For 
a continuum of levels, such as in metallic quantum dots, 
A = 0, incoherent cotunneling dominates and we expect 
little asymmetry. Dephasing due to electron-electron in- 
teractions in the dot is-beyond the scope of the present 
analysis, but is knowncj to be inefficient. 



(D 






FIG. 6. Sequence of asymmetry of AB oscillations in the 
Coulomb blockade regime. The solid line depicts schemati- 
cally the conductance oscillations vs. the gate voltage V g , and 
shows the Coulomb peaks. Regions of small and large AB am- 
plitudes are indicated, showing the asymmetry between two 
adjacent Coulomb blockade valleys ("large" corresponds to a 
valley with a total spin 0, while "small" corresponds to an 
S = 1/2 valley). 

To investigate how the crossover occurs we count the 
number of coherent and incoherent cotunneling channels 
for a system with arbitrary level spacing A. We examine 
four different situations. In all of them the level energy e 
is tuned away from resonance, AsT, T <C |e| < U/2. The 
four situations arise due to the combination of being in 
a valley which corresponds to an even or odd electron 
number of the dot and being on either side of the res- 
onance, i.e., only particle- like or hole- like cotunneling is 
considered. 

When the level spacing A becomes comparable to |e| 
the number of dot levels involved in transport is increased 
from 1 to N ~ 2|e|/A + 1 (each of them is spin degen- 
erate). As long a£ A > fc^T, no inelastic cotunneling 
process takes place which changes the energy of the QD. 
Incoherence occurs only due to spin flip. In the "even" 
valley all channels are coherent while in the "odd" valley 
1 out of N + 1 corresponds to spin flip and is, there- 
fore, incoherent (the total number of channels is always 
the same) . The relative asymmetry around the resonance 
peaks is N/(N+1) and vanishes for large N, i.e., A <C \e\. 

Now we reduce A further such that A <C fc^T. The 
number of levels within the range defined by temperature 
is M ~ 2/cbT/A ^> 1 (each of them is spin degenerate). 
For a simple estimate we assume that all single-particle 
states within the energy range ifc^T around the Fermi 
energy in the leads are occupied with probability 1/2, all 
states above that range are empty, and all states below 
arc filled. For the single-particle states in the QD we do 
the same. Let S be the level spacing in the leads with 
5 -C A. The number of coherent channels (times the 
probabilities that the corresponding states are empty or 
occupied) is then 



A NM 



(3.20) 



spin of the incoming and outgoing electron there are 
[M + 3(E/A)]/4 combinations to leave this excitation 
on the QD (where the probabilities to find the corre- 
sponding state empty or filled are already included). For 
E < (i.e., energy is pulled out of the QD) we have 
[M + (E/A)]/4 combinations. Energy conservation re- 
quires that the energy stored in the QD is pulled out 
of the leads and vice versa. An analogous considera- 
tion yields for E > (energy pulled out of the leads) 
(A/S)[M + (E/ A)]/4 combinations and for £<0we get 
(A/8)[M + 3(E/A)]/4. In total we find, therefore, 

5 • T ' X> + 3n ^ M + ") = 4 [ Af3 + ° ( Af2 )] 



2 S 



(3.21) 



Next, we count the number of incoherent channels 
where an excitation of energy E is left on the QD (E 
is an integer multiple of A). For E > and given the 



channels with either E > or E < (a multiplication 
factor of 4 was introduced to account for the spin of the 
incoming and outgoing electrons). The number of inco- 
herent channels with E = (and spin flip) is of the order 
(A/S)M 2 and can be neglected. 

The ratio of the number of coherent channels (for 
N ^> M ^> 1) to the total number of transmission 
channels is N/(N + AM 2 ). As a consequence the co- 
herent contribution vanishes for ksT ^ -y/|e| A./8 or 
8(fc B T) 2 /|e| » A. 

In summary, by reducing the dot level spacing A we 
go through three different regimes: 

(i) A> |e|: the asymmetry due to spin-flip cotunneling 
is as large as possible; 

(ii) |e| ^> A ^> 8(fcsT) 2 /|e|: the asymmetry vanishes but 
transport is dominated by coherent channels; 
(iii) 8(fcsT) 2 /|e| ^> A: transport is dominated by inco- 
herent channels. 



F. First-order transport 

The discussion so far covered cotunneling, which dom- 
inates away from resonance. Let us now turn to the 
regime fc^T ^S> I\ |e|, where the dot level is near reso- 
nance and transmission is dominated by first-order trans- 
port in r. Can we tune the single-dot AB interferometer 
in such a way that full destructive interference is achieved 
for all those first-order transport channels that are coher- 
ent? The energy spread of electrons going through the 
reference arm is fcgT, while the width of the resonance 
through the QD is V. Hence, the temperature has to be 
on the one hand larger than T in order to be in the regime 
where first-order transport dominates, yet, on the other 
hand, it has to be smaller than T to allow for a destruc- 
tive interference of all energy components simultaneously. 
To circumvent this problem, we consider a two-terminal 
AB interferometer with two QDs, one in each arm, see 
Fig. Be. 



10 



IV. INTERFEROMETRY WITH TWO QDS 

The conceptual difficulty to address first-order trans- 
port in a single-dot AB interferometer is that the temper- 
ature has to be on the one hand large, yet, on the other 
hand, it has to be small to allow for a destructive interfer- 
ence of all energy components simultaneously. As we will 
show in this section, this difficulty will not arise in AB 
interferometers which contain two QDs, one in each arm, 
see Fig. 0c. This is due to the fact that the resonance 
of width r of each QD filters out a small fraction of the 
incoming electrons in both arms, even at high tempera- 
ture, and we will find that fully destructive interference 
(in the absence of interaction) is feasible. For interact- 
ing QDs we find again an asymmetry factor which makes 
it possible to distinguish on which side of the resonance 
the adjacent Coulomb blockade valley has an even or odd 
number of electrons. 

Geometries similar to the two-dot AB interferome- 
ter have been investigated theoretically in the literature. 
Resonant tunneling (in the absence of interaction)^ and 
cotunnelingJdO have been studied in this geometry. Fur- 
thermore, a numerical renormalization-group analysis 
has been performed!! where, by construction, only equi- 
librium properties of the dot were included. Spectral 
properties of such a double-dot systeiBp-pr equivalently 
a two-level dot, have been addressed. E3e3 Recently, also 
experimental-realization of a two-dot AB setup has been 
reportedMN 



A. General current formula 



The total Hamiltonian 



H = H 



dot,l 



rrdot,2 



(4.1) 



is the sum of two parts.c3 Each of them describe a QD 
coupled to (the same) reservoirs and has the structure of 
Eq. (2.1). We choose a completely symmetric geometry, 
and we assume ksT ~^> V, |ei|, |e2 1 as well as T ^ |ei — ea|, 
where ei,2 is the energy of the level in QD 1 and 2. In this 
regime lowest-order transport dominates, and we can set 
e = 6\ = e 2 . To model the enclosed flux we attach a phase 
factor e iv ' 4 to the tunnel matrix elements i^doti and 
ti,dot2, and e~ lip/i to tL,doti and £fl,dot2- This symmetric 
flux dependence of the tunnel matrix elements can be 
achieved choosing a corresponding gauge. 

The system is equivalent to a single QD with two levels 
(each of them spin degenerate) with (^-dependent tunnel 
matrix elements. To write the total current in a compact 
way we employ &J2~x 2 matrix notation to account for the 
two QDs and getH 

I tot = ^Jdutr{[T L f L -T R f R }G> 

+ [T L (l-f L )-T R (l-f R )]G<} (4.2) 



1 



o+if/2 



^h I* = ^ e _t v/2 e x J <W and I* = (I*) . 

The trace has to be performed over both the spin degrees 
of freedom and the 2x2 matrices. We get the general 
and exact result 



Oh /-* 



du 



2ft 

+ cos|(Gf 2>CT 



22, a- "22,(7 



"12,<7 



(h - Ir) 

G 21,<7 - G 21,a) (/i ~ I r) 



-2i sin ■ 



f ( Il + Ih 



-2i sin 



(""12,(7 "21,o 



1- 



h + Jr 



(^12,(7 — ^21,, 



(4.3) 



with the full Green's functions which includes multiple 
tunneling through the entire device. The indices 1 and 
2 label the respective QD. Note that even the diagonal 
terms Gn ]dr and G 22>a are not the Green's functions of 
QD 1 and 2 in the absence of the other one. Instead, 
they have to be calculated in the presence of the entire 
geometry. 



B. First-order transport in linear response 

To address the low est-order contribution to transport 
we expand Eq. ( |4.3| ) up to first order in T, i.e., only 
Green's functions in zeroth order are involved. The off- 



>(o), 



-»<(0)/ 



diagonal terms are connected by G 12 : o :(w) = G 12 J(u) — 
2iriP 2 £S(u> — e) to the stationary off-diagonal density- 
matrix elements P\° = (\2a)(la\) in zeroth order in 
r. In equilibrium (and zeroth order in T) the density 
matrix is diagonal, with the probabilities determined 
by Boltzmann weights, and all off-diagonal matrix ele- 
ments vanish. In the first and second line of Eq. flf.3| ), 
only equilibrium Green's function are involved in linear- 
response transport. As a consequence, the second line of 
Eq. (4.3) does not contribute at all (in linear response 
and first order in T). In the third and fourth line, how- 
ever, nonquilibrium Green's functions enter even in the 
linear-response regime. We obtain 



OP 



dV 



= ^rW^{lmG^M/'(, 



.■ <P 



dG 



>(0) 



+ sin — f(u>) ■ „ 

2 Jy ' d(eV) 

a r <(0) 

+ Sm 2 [1 - /H] 9P0 



(4.4) 



Here, wc have used the fact that the contributions in- 



-,r(0) 



,<(0) 



>(0), 



volving ImG 2 2j(u), G 2lCT (w), and G 2 iJ(uj) amount 



to an overall factor 2. For the first term in Eq. (4.4) 



-,r(0) 



we use — (l/7r)ImG 1 ^ a {uj) = 8{uj — e) for U — and 
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-(1/tt) ImG^liuj) = %-e)/[l+/(e)] for U = oo. This 
term is flux independent and is twice the conductance 
through one QD in the absence of the other one. Inter- 
ference effects are accounted for by the second and third 
terms. To determine the off-diagonal density-matrix el- 
ement P\° we use the real-time transport theory devel- 
oped in Refs. |2^ and |34| and solve a generalized master 
equation. The analysis, presented in Appendix pi em- 
ploys diagrams of the type used in Appendix |A|. We find 
the result (for zeroth order in T and V = 0) 



3{eV) 



-/'(e) sin^/2) 



in the absence of interaction and 



dP 



• la 



/'(e) 



d(eV) ~ 2 [1 + /(e)]* 



i(tp/2) 



(4.5) 



(4.6) 



for U = oo. As a consequence, in the absence of a n AB 
flux, only equilibrium Green's functions enter Eq. (4.4). 
It is crucial, however, that in the presence of the flux, 
nonequilibrium properties of the Green's functions are 
involved even for the linear-response regime. 

Collecting all terms we find for the noninteracting case 



dP 



dV 



or 



dot 



V=0 



dV 



v=o 



[1 



W2)] (4.7) 



with (dI dot /dV)\ v =a = -{ire 2 /h)rf(e) being the con- 
ductance through a single QD. At tp = ±7r, ±37r, . . ., the 
total current vanishes (up to a small correction of or- 
der |ei — £2|/r, which is not considered here), indicating 
that lowest-ordex-.transport is fully coherent. In fact it 
has been shownliJ that in this case, the system can be 
mapped onto a model with two levels, where each level 
is coupled to one of the leads only. As a consequence, in 
the absence of any interaction, there is no way to transfer 
any electron from one lead to the other. This completely 
contrasts the picture of sequential tunneling which pre- 
dicts incoherent transmission. In the absence of inter- 
action, however, the transport should be fully cohe rent . 
For the simple limit U = 0, we can rederive Eq. (4.7) 
by using the Landauer-Biittiker approach, determining 
the dot Green's function by equations of motion, and ex- 
panding the result up to first order in I\ This is done in 
Appendix |(J. The sequential-tunneling picture is highly 
misleading in this context. We should rather view the 
transport as resonant tunneling but expanded up to low- 
est order in I\ 

In the presence of interaction we obtain 



dP 



dV 



= 2- 



81 



dot 



V=0 



dV 



v=o 



1- 



sin 2 (<^/2) 



[l + /(e)F 



(4.8) 



with (dI d ° t /dV)\ v =o = -(ne 2 /h)Tf(e)/[l + /(e)]. 

We point out that the maximal total conductance 
(which is reached at tp = 0, ±2ir, ±47r, . . .) is the sum of 
the conductances through the QDs taken apart. There is 



no extra factor 2 familiar from constructive interference, 
which takes place away from resonance when transmis- 
sion is small. To understand this we remark that for 
tp = only the symmetric combination of the dot levels 
are coupled to the leads while the antisymmetric one is 
not (this holds for degenerate levels which are symmet- 
rically coupled to the leads). At U — the model is, 
therefore, equivalent to a system with one level but tun- 
nel matrix elements which are y/2 times the original ones. 
As a consequence, in first order in T (which is second or- 
der in the tunnel matrix elements) the total conductance 
is twice the conductance through one arm of the inter- 
ferometer in the absence of the other one. 

The factor 1/[1 + /(e)] 2 yields an interaction-induced 
asymmetry in the ratio of coherent to total transport 
around a conductance peak. This factor can be easily 
understood in the following way. Consider the situation 
where the incoming electron has spin up. Interference 
is possible if we start with both QDs being empty, one 
QD empty and the other one filled with spin up, or both 
QDs filled with spin up. For all other starting configura- 
tions one would be able to tell afterwards which way the 
transfered electron has taken. The corresponding con- 
tributions are not subject to interference. Summing up 
the probabilities for the four situations which allow for 
interference we find (at V = 0) the factor 1/[1 + /(e)] 2 . 

The effect of the interaction-induced asymmetry factor 
is illustrated in Fig. 0. 




FIG. 7. Asymmetry in interference signal for AB interfer- 
ometer with two QDs. Main panel: maximal conductance 
(i.e., tp — 0) as a function of level position. Inset: AB os- 
cillations for the two level positions indicated by a and b in 
the main panel. All conductances are normalized by F /(fcsT) 
and plotted in units of e 2 /h. 

In the main panel the maximal first-order linear con- 
ductance (which occurs at tp = 0) is displayed. In our 
simple model the resonance peak is symmetric around 
e/iksT) = (1/2) In 2. Now we pick two level positions 
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at which this maximal conductance is equal, indicated in 
the figure by a and b. The amplitudes of the AB oscilla- 
tions at these points differ strongly from each other (in- 
set of Fig. 0), due to the asymmetry factor 1/[1 + /(e)] 2 . 
The smaller AB amplitude corresponds to the situation 
in which single occupation of each QD is more likely than 
the dots being empty. 



C. Effect of dot-dot interaction 

We want to emphasize that for a proper realization of 
an AB interferometer which might provide full destruc- 
tive interference of coherent transport channels, the two 
QDs should not influence each other electrostatically. Let 
us, e.g., consider the case of spinless electrons in an AB 
interferometer in which simultaneous occupation of both 
QDs costs the elect rostatic-energy penalty U'. Again we 
start with Eq. (4.4), solve for the corresponding density- 
matrix elements (see Appendix pf), and end up with 



dP 



dV 



dl 



dot. 



V=0 



dV 



x [l-Csin 2 (vV2)] (4.9) 



v=o 



with 



C = 



1 



\n 2 ((3U'/2ir)cos 2 (<p/2) 



! [l"/(e 



(4.10) 



for 0e ~ 1 and f3U' > 1. At tp = 0,±27r,±47r, . . ., the 
total linear conductance is the sum of the conductances 
through the dots taken apart. For tp = ±7r, ±37r, . . ., the 
conductance is zero. The latter statement is not only 
true for first-order transport but a general fact, since it 
is possible to rotate the basis of the double-dot system 
such that one state only couples to the left reservoir and 
the other one to the right reservoir .E3 

We observe that now, since the two paths from source 
to drain are influencing each other electrostatically, the 
interference signal no longer provides a direct tool to dis- 
tinguish coherent from incoherent transport through a 
single QD. The amplitude of the oscillations is reduced 
for increasing U' because an electron occupying one QD 
effectively blocks the path through the other QD, and 
interference is suppressed. 



V. FINITE-BIAS-INDUCED AB OSCILLATIONS 



described in zeroth order.) Out of equilibrium, however, 
a flux dependence of (N) might arise, as we will show in 
this section. 



A. Single-dot AB interferometer 

In the absence of charging energy, the correction terms 
of the probabilities for the empty (x = 0), singly- 
occupied (x — & with a =] or {), and doubly-occupied 
dot (x = d), in first order in i rcf are (see Appendix ^) 



P { x ] = a xY^^\ fCf I sin ^ IMe) ~ Me)] 



(5.1) 



with oo = 1 - F(e), a a = F(e) - 1/2, and a d = -F(e). 
Here we have used the definition 



F(e) 



T L fL(e)+T R f R (e) 

Tr+Tn. 



(5.2) 



At finite bias voltage, the probabilities for the dot states 
depend on the AB flux even in zeroth order in the in- 
trinsic linewidth T. Only at special values of the flux, 
tp = 0, ±7T, ±2tt, . . ., the correction terms P^ vanish for 
all states \ — 0, T> -L d, and the probabilities coincide with 
those of a QD without a reference arm. 

The same statement is true for interacting QDs. For 
U = oo we obtain correction terms of the type Eq. (5.1) 
with a® = 1/[1 + F(e)] 2 and a a = —o.q/2 (see Ap- 
pendix ^) . 

For symmetric couplings, Tl = Fr, the occupation of 
the QD up to first order in t rcl is 

(N) = f L (e) + f R (e) - \f cf \[f L (e) - / R (e)] sin^ (5.3) 
for noninteracting QDs, U — 0, and 



(N) = 



U(e)+Me) |i ref |[/L(e)-Me)]sin^ 



1 + 



MibfeH 



1 



fL(e)+f R (e) 
1 



(5.4) 



for U = oo. 

The generation of AB oscillations due to finite bias is 
illustrated in Fig. ||. To be specific we chose the level 
energy e where the conductance has its peak, e = for 
U = and e/(k B T) = (1/2) In 2. Furthermore, we used 
\t Icl | = 0.1 and chose Tl — Fr in this example. 



In this section we concentrate on the average number 
(N) of electrons in the QD near resonance, |e| < ksT, in 
the weak-coupling regime, T <g; ksT (and r <C A). At 
equilibrium, (N) is, for both the single-dot and two-dot 
AB interferometer, determined by classical Boltzmann 
weights. The latter depend on the energy of the dot level 
only. In particular, they are independent of magnetic 
flux. (In the regime specified above, transport is domi- 
nated by first order in T, while the occupation number is 
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0.4 0.6 0.8 1.0 

(p/27i 

FIG. 8. Bias-voltage-induced AB oscillations of the dot oc- 
cupation in a single-dot AB interferometer. The level position 
is tuned at the maximum of the Coulomb-oscillation peak. In 
equilibrium the average charge on the QD is flux independent 
(in lowest order in V). With finite bias voltage, AB oscil- 



lations emerge. 

r L = r«. 



In this example we chose |i rc | = 0.1 and 



B. Two-dot AB interferometer 

For the AB interferometer containing two QDs it turns 
out that the corrections to the diagonal matrix elements 
of the density matrix in lowest order in T remain zero 
in first order in V, i.e., in contrast to the single-dot AB 
interferometer, there is no bias-voltage-induced AB os- 
cillation of the occupation of either dot in linear order in 
V. This changes when electrostatic inter-dot interaction 
is introduced (although it does not change the spatial 
symmetry). For the model with spinless electrons and 
inter-dot interaction we find that the correction for the 
occupation of QD 1, P*, is determined (for /3e ~ 1 and 
f3U' > 1) by 



dPl 



dV 



CH(3U'/2tt) /'(e) 



v=o 



47r[l-/(e)P l + /(e) 



sm ip , 



(5.5) 



where t he ( also flux-dependent) term C is defined 
in Eq. (p~To|) . For QD 2 we get (dPi/dV)\ v=0 = 
— (dPl /dV)\v=o- All these results are derived in Ap- 
pendix H. 



VI. BREAKDOWN OF PHASE LOCKING AND 
MORE GENERAL SYMMETRY RELATIONS 

Onsager relations yield phase locking of the linear con- 
ductance through a (closed) two-terminal device, i.e., the 
linear conductance is symmetric under reversal of mag- 
netic flux, p — > —ip. Phase locking is, however, no longer 



enforced in the nonlinear-response regime. In the final re- 
sult for the interference current through a sing le-dot AB 
interferometer in first order in \t \, Eq. (3.12), the first 
line is symmetric under ip — > —ip, but the (^-dependence 
of G^ in the second line can break this symmetry at fi- 
nite bias voltages. In the following we explicitly calculate 
the current in lowest order in T. It will turn out that for 
a noninteracting QD phase locking survives even at finite 
bias, but is broken for interacting QDs. 

In the noninteracting case we find to lowest order in T 



that ImG. 
yields 
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/;." = -rv / fIf^|i r 



| COS If I did 



/l(w) -/r(w) 



(6.1) 



where the prime at the integral indicates Cauchy's prin- 
cipal value. This proves that for a noninteracting QD 
phase locking is preserved even in nonequilibrium, at 
least in lowest order in T. It can be shown, however, 
that even in higher orders in T the Green's function is 
invariant under ip «-> —ip, i.e., phase locking is preserved 
in nonequilibrium. This is due to the fact that the to- 
tal current is obtained as a sum over contributions at 
different energies, with weights given by the nonequilib- 
rium conditions. Each of these contributions satisfies the 
phase-locking symmetry separately. 

This changes when interaction in turned on. For U = 



oo we find that Im G 
reads 



r(l) 



is non zero, and the current 



r(l)_ 
1 r? — 



4eyTZr^|* rcf |coBp 
h 1 + F(e) 

, 4^e (T L T R ) 3 / 2 \f 



r' h(w) 

auj 



/*(w) 



| simp 



r 2 [l + F(e)] 2 



lh(e)-f R (e)] 2 . (6.2) 



In addition to the first line, which corresponds to phase 
locking, there is a contribution proportional to sin ip (sec- 
ond line) which breaks the symmetry under ip <-> —ip. As 
we see it takes both interaction and finite bias voltage to 
break phase locking. 

For the AB interferometer which involves two QDs, 
phase locking remains at finite bias, independent of 
whether the interaction is included or not. This is due to 
the symmetric setup we chose and follows from general 
symmetry relations, as discussed next. 



Connection between spatial symmetry and 
symmetry of transport coefficients 

The general relation for all two-terminal setups 
dI(V,<p) dI(-V,-<p) 



dV 



dV 



(6.3) 



where V is the symmetrically applied bias, yields as a 
direct consequence the Onsager relation 
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dlfr) 



dV 



di(-cp) 



v=o 



dV 



(6.4) 



ai(y,<p) di(-v,-<p) 



v=o 



which requires phase locking in linear resp onse. 

In addition to this symmetry Eq. ( |6.3| ), spatial sym- 
metries can be exploited as well. Let us first consider an 
AB interferometer with a single QD. If T^ — Tr then the 
system remains invariant if it is mirrored with respect to 
a vertical axis and the chemical potentials of the left and 
right leads are exchanged. This leads only to Eq. ( |6.3| ) 
and does not induce a new relation. 

The situation is more complex for the AB interferom- 
eter made out of two QDs. We assume that the QDs are 
identical. The coupling of the QDs to the leads, how- 
ever, may be different. We consider three different cases 
in which the system has a spatial symmetry, see Fig. ^[ 



a) 



QD1 



H+eV/2 



H-eV/2 



b) 



c) 




|U+eV/2 



(i-eV/2 



H+eV/2 



H-eV/2 



QD2 

FIG. 9. Two-dot AB interferometers with different spatial 
symmetry. Thick and thin lines indicate stronger and weaker 
tunnel coupling, respectively, a) possesses mirror symmetry 
with respect to a vertical axis, b) is invariant under rotation 
at angle it, and c) has mirror symmetry with respect to a 
horizontal axis. In cases b) and c) phase locking is preserved 
even for finite bias. 

If kL,doti| = |*R,doti| and |tx,,dot2| = |<R,dot2| (see 
Fig. KJa) the system has mirror symmetry with respect 
to a vertical axis only, and the resulting symmetry rela- 
tion 



dV 



dV 



(6.5) 



coincides with Eq fl6.3p . 

If liL,doti| = |tfl,dot2| and |t_R,doti| = |*L,dot2| (see 
Fig. |9p) the system is invariant under rotation at angle 
■k in the plane and reversal of V . As a consequence 



di(y,<p) di(-v^) 



dV 



dV 



(6.6) 



If J*L,doti| = |^L : dot2| and |tfl,doti| = |*«,dot2| (see 
Fig. Bp) the system has mirror symmetry with respect 
to a horizontal axis. This implies 



ai(y,(p) di(v,-<p) 



dV 



dV 



(6.7) 



In the two latter cases phase locking occurs, wh ich either 
follows directly or after making use of Eq. (6\3). This is 
a consequence of spatial symmetry. In the first case, or 
without any spatial symmetry, breaking of phase locking 
may be possible for finite voltages. 



VII. SUMMARY AND DISCUSSION 

The focus of this work is the interplay between coher- 
ent and incoherent transport channels of a mesoscopic 
setup. We specifically address the issue of transport 
through a quantum dot. To stress the difference be- 
tween dephasing processes (which do not require any en- 
ergy exchange) and inelastic scattering (which, of course, 
may also introduce incoherent transport), we have cho- 
sen to discuss spin-flip processes as the source of dephas- 
ing. Our analysis then elucidates some important facets 
of electronic spin transport. Central to our analysis is 
the fact that the existence of decohering spin-flip chan- 
nels requires the presence of electron-electron interaction 
(in our model the constant capacitance term). We have 
shown explicitly that in the absence of such interaction 
the contributions from the spin- flip channels cancel out. 

In order to quantify and study the "degree of coher- 
ence" in the system we had to resort to interferomctry 
setups. This is why we have embedded our system in an 
Aharonov-Bohm circuit, in line with recent experiments. 
As we present results for the AB amplitude of the trans- 
mitted current (as compared with the flux-insensitive 
component of the current), our results are closely related 
the "visibility" , employed by some experimentalists in 
this context. 

Technically we have performed a systematic analysis 
of first- and second-order (in the QD-lead coupling, T) 
transport through AB interferometers containing cither 
one or two QDs. As a consequence, in the Coulomb 
blockade valley with one dot level being singly-occupied 
only one half of the low-temperature (yet T > Tk) trans- 
port is coherent. This is to be contrasted with the case 
of all dot levels being either empty or doubly occupied, 
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where transport is fully coherent. The conclusions of 
our analysis, performed for a QD with a single level, are 
generalized to a multilevel QD. In contrast to expecta- 
tions based on a sequential-tunneling picture, we have 
proved that even in first-order transport, the transmis- 
sion is at least partially coherent. Our formalism allows 
to cover the nonlinear-response regime, too. We predict 
bias-voltage-induced AB oscillations in the occupation 
number of the QD of a single-dot AB interferometer. 

An important outcome of our analysis is the partial 
suppression of the AB conductance oscillations which, in 
the presence of e-e interactions, is asymmetric with re- 
spect to the position of the Coulomb peaks. This asym- 
metry is a manifestation of the presence of incoherent 
transport channels, whose number is different between 
the two sides of a Coulomb peak. Our asymmetry effect 
(calculated in low orders in T) increases with lowering 
the temperature. The physical picture is expected to 
change, though, as the temperature is lowered towards 
and beyond the Kondo temperature Tk- In this case, 
the local spin of the QD is screened by the lead electron 
spins. At zero temperature the QD is seen by the trans- 
mitted electrons as a pure static potential scatterer, i.e., 
no spin-flip appears. We therefore expect an increase of 
the coherent transport contribution when the tempera- 
ture is decreased below Tk, and a gradual disappearance 
of our asymmetry effect. This shows that while our co- 
herence asymmetry effect may provide information com- 
parable to that obtained through the Kondo physics (i.e., 
the QD's spin), the latter reflects different physics. We 
also note that the asymmetry we refer to concerns the 
magnitude of the AB amplitudes at two points on dif- 
ferent sides of the Coulomb blockade peak, which cor- 
respond to the same total (flux-averaged) conductance. 
Our asymmetry does not imply an asymmetric line shape 
of the conductance, and therefore is different from the 
Fano mechanism. The presence of such asymmetries in 
the experiment will be discussed elsewhere.H The results 
of our analysis suggest that one may use this asymme- 
try to probe the total spin of a QD; we discuss how this 
asymmetry is suppressed for multilevel QDs, as the level 
spacing is reduced. 

We have included in our analysis a brief discussion of 
time-reversal symmetry in the problem, specifically the 
issue of phase locking. The latter is usually presented as 
an outcome of Onsager relations. We have shown how, in 
general, it breaks down at finite bias voltages. In certain 
cases, however, we find that phase locking is preserved 
due to spatial symmetry of the setup. 
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APPENDIX A: DENSITY MATRIX FOR 
SINGLE-DOT AB INTERFEROMETER 

In this appendix we present the somewhat technical 
part of calculating the density-matrix elements for the 
single-dot AB interferometer. The derivation is based on 
the real-time transport theory developed in Refs. E9 34. 
The starting point is the generalized stationary master 
equation in Liouville space, 



(e> 



pXi 

Xl 



E 

X\_ --X-2 



pXiyOCi,Xt 

X2 X2'X2 



0. 



(Al) 



where \i an d X2 denote any state for the QD, e Xl and 
e X2 are the corresponding energies, and P* 1 = (|x2)(xi|) 
is a matrix element of the reduced density matrix for the 
QD subsystem. The diagonal matrix elements P Xl = P* 1 
are nothing but the probabilities to find the system in a 
given state xi- 

The matrix elements P of the density matrix are con- 



nected to each other in Eq. (Al) by the terms £ 



x 2 .X2 

which can be viewed as generalized transition rates in 
Liouville space. They are defined as the irreducible self- 
energy parts of the propagation in Liouville space and 
are represented as diagram blocks on a Keldysh contour. 
For a detailed derivation of this diagrammatic language, 
the generalized master equation, and the rules on how to 
calculate the value of a diagram we refer to Refs. 29 33. 

In the present context we calculate the corrections lin- 
ear in t to the probabilities for the dot states. Wc 
concentrate on the regime where the dot level is tuned 
close to resonance, i.e., we determine the probabilities in 
zcroth order in T. To this end we need the self-energy 
parts E to lowest (first) order in T. The possible dot 
states x are labeled by for an empty QD, a for single 
occupancy with spin a =|, J,, and al for double occupancy. 

First, we o bserve that in the generalized master equa- 
tion Eq. (Al) the probabilities P x = P* couple to diag- 
onal matrix elements only. This has to do with the fact 
that the z component of the spin is a conserved quan- 
tity. As a consequence, only self-energy parts of the type 
S XiX ' = "£> xx i enter. We get to zeroth order (indicated 
by the superscript (0)) in t ref the master equations 
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together with the normalization condition Pq ' + 2P„ ' + 

P d = 1. The diagrammatic representation of the irre- 
ducible self-energy parts are shown in Fig. HO. 
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FIG. 10. Irreducible self-energy parts to zeroth order in t Te 
(and first order in V). 

There are two horizontal lines representing the forward 
(upper line) and backward (lower line) propagators on 
the Keldysh contour. Vertices (full dots) indicate tun- 
neling from the QD to the lead or vice versa. They are 
connected in pairs by tunnel lines (dashed lines) repre- 
senting contractions of lead electron operators. Depend- 
ing on the order of the annihilation and creation opera- 
tors, specified by the arrows, the tunnel lines contribute 
with a Fermi function f r (if the line goes backward with 
respect to the Keldysh contour) or 1 — /,. (if the line goes 
forward), r = L/R, times |t r | 2 iV r = r. r /(27r), where N r 
is the density of states in reservoir r. Along the tunnel 
lines and the forward and backward propagators we as- 
sign the proper dot states such that the particle number 
and z component of the spin are conserved at each vertex. 

It is easy to see that it is impossible to construct a 
diagram for Udo or Eod to first order in T, i.e., E^o = 
Eod = 0. The other diagrams are depicted in Fig. [f(J. 
As an example we calculate E^q explicitly using the di- 
agrammatic rules derived in Refs. p9,£54l The tunnel line 
contributes with rV/(27r)[l — / r (w)]. There is a minus 
sign for each vertex on the backward propagator. This 
yields a total minus sign for the diagram under consider- 
ation. For each segment between two vertices there is a 
resolvent 1/(AE + i0 + ), where AE is the difference be- 
tween the energy on the left-going lines and right-going 
lines (including the propagators and the tunnel lines). 
For the left diagram we obtain l/(e — u> + i0 + ). Finally, 
we integrate over the energy u) of the tunnel line and sum 
over the reservoir index r to get 
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The value of the right diagram is the same except for the 
denominator which now reads u> — e + i0 + . After making 
use of l/(x + i0 + ) = P(\/x) — iwd(x) we can evaluate 
the integrals. The final result (for all diagrams shown in 
Fig. M is 
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where we have used the definitions F(uj) = [rx,/i(w) + 
T R f R (uj)}/(r L + T R ) and r = T L + T R . This yields the 
solution 



P 



(„)_ [1-F(e)][l-F(e + U)} 

F(e) + 1-F(e + U) 
(0) = F{e)[l-F(e + U)] 
F(e) + 1-F(e + U) 
(0) = F(e)F(e + U) 
d F(e) + 1-F(e + U)' 



(A9) 
(A10) 
(All) 



In equilibrium P(u>) = f(co), and we recover the clas- 
sical Boltzmann factors which determine the probabil- 
ities, P (0) = 1/Z, p] 0) = exp(-/3e)/Z, and P^ 0) = 
exp[-0(2e+U)]/ZwithZ = l+2exp(-/3e)+exp[-/3(2e+ 
[/)]. 

The effect of the reference arm on the QD shows up in 
the correction terms to the probabilities in first order in 
t lci . The master equations to solve are 

^ocV + 2Pi°^ + P„ (1) 4° > + 2P«4°> = (A12) 
P (o) E (i) + 2P (0) E (i) + P (D S (0) + 2p ( 1)E (0) = (A13) 

along with the normalization condition P + 2P^ + 
0. The diagrams for the self-energy part to first 
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order in t rcf are shown in Fig. [Ll| 
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FIG. 11. Irreducible self-energy parts to first order in t Tci 
(and first order in V). There are two kinds of vertices: single 
vertices (full dots), which represent tunneling between QD 
and leads, and double vertices (open dots), which indicate 
direct tunneling between the leads. We show explicitly all 
diagrams with direct tunneling from left to right (they are 
proportional to e' v ). The remaining diagrams are obtained 
by exchanging the labels L and 7? for the left and right lead 
everywhere (those are proportional to e~ 1 ^). 
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Here, another type of vertex has to be introduced 
which represents direct tunneling from one lead to the 
other. Such a vertex is connected to two tunnel lines 
("double vertex") since it represents two lead electron 
operators which are contracted with another operator 
of the same lead each. As a consequence we encounter 
objects made of two tunnel lines connecting one double 
vertex with two single vertices. Each line contributes 
with f r or 1 — f r , and the total multiplication factor 
is N L N R t L t R \i\ = |i rof |Vrro(4^ 2 ) times e±^. The 
upper (lower) sign applies if the incoming line into the 
double vertex belongs to the left (right) lead. 

Again we evaluate one diagram explicitly, namely the 



.(i) 



first one contributing to E^g as shown in Fig. [11 
find 



We 



duiL / dtoft 



(wl - Wr + iO+)(e -C0 R + i0+) 



(A14) 



times \f ci \e tv yJT lT R / (Ait 2 ) . There is one minus sign due 
to the vertex on the lower propagator. This is, how- 
ever, canceled by another minus sign coming from the 
ordering of the Fermi operators: in the diagrammatic 
language, each crossing of tunnel lines introduces a mi- 
nus sign (to apply this rule view the double vertex as two 
separate vertices lying close together, where the vertex 
connected to the incomiaff line comes first with respect 
to the Keldysh contour). EJ 

We note that for all self-energy parts shown in Fig. [ll] 
the fifth and sixth diagrams cancel each other out since 
they differ by a minus sign due to the vertices on the 
lower propagator. The same holds for the eleventh and 
twelfth diagrams. 

After collecting all contributions we end up with 

E« = i sin ¥ VrIr^|i ref | [f L (e) - f R (e)] (A15) 

4d = ^sm^T L T R \f cf \ [f L (e + U)~ f R (e + U)} (A16) 

as well as E$ = 2E« and E« = 2E« 

The results for P x ' are lengthy expressions. They sim- 
plify, however, for either U — or U — oo. For these 
limits they are given by Eq. (|5.l|). 



APPENDIX B: DENSITY MATRIX FOR 
TWO-DOT AB INTERFEROMETER 

In this appendix we determine the off-diagonal matrix 
element P 2 1 <t m zeroth order in T for the two-dot AB in- 
terferometer. These results are needed for the evaluation 
of Eq. (| 



To achieve this goal we expand Eq. (Al) up 
to zeroth and first order in T. The irreducible self-energy 
parts E have contributions of order T and higher. As 



a consequence, the zeroth order expansion of Eq. (Al) 
yields that all off-diagonal matrix elements P* 1 with 
Xi 7^ X2 c a n only arise for e Xl 
vanishes. 



The calculation in this section is based on the same 
diagrammatic language as introduced in the previous ap- 
pendix. In the two-dot AB interferometer there is no di- 
rect tunneling between leads, i.e., no double vertices will 
appear in the diagrams. Since constructing and evaluat- 
ing the diagrams is straightforward along the line indi- 
cated for the single-dot AB interferometer in the previous 
appendix, we refrain from drawing them explicitly. 



1. No interaction 

In the case of noninteracting QDs, [7 = 0, the two 
channels provided by the spin degree of freedom do not 
influence each other. The Hamiltonian is just the sum of 
two identical models, one for each spin. In this case it is 
sufficient to consider a simpler model with spinless elec- 
trons and multiply the final expressions for the current by 
a trivial factor of 2. This simplifies both the calculations 
and the notations. 

The Hilbert space of the double-dot system is, then, 
four dimensional: the double dot can be empty (% = 0), 
singly occupied with the electron in dot 1 or 2 (x = 1, 2), 
or both dots are filled (x = 12). The corresponding ener- 
gies are 0, e, and 2e, respectively. The density matrix has 
16 matrix elements. Four of them are the diagonal ma- 
trix elements Pq , Pi , P 2 , and Pjf • They are all real and 
positive. Since the transformation 1 <-» 2 and V *-> —V 
does not change the system, the three combinations P§, 
Pi + P 2 2 , and Pjf are even in V, while the combination 
Pi — P 2 2 is odd. There are two nonvanishing off-diagonal 
matrix elements, P 2 and P2 1 . All others are zero in zeroth 
order of T since the energies of the corresponding states 
differ from each other. The hermiticity of the density 
matrix implies P 2 = (P2 1 ) ■ 

In equilibrium, V = 0, the diagonal matrix elements 
P* are probabilities for the state x determined by clas- 
sical Boltzmann factors exp(— /3e x ), and all off-diagonal 
matrix elements vanish. Our goal is to determine the lin- 
ear corrections in V. Only the elements which are odd in 
V have a finite correction. These elements are P/ — P$ 
and Im P 2 . 

In the following, we write P = P + P + . . . and E = 
E + E + . . . for the zeroth and first order in V. It turns out 
that there are two independent equations which relate the 
first-order corrections Pi = — Pf and ImP 2 1 = — ImP 2 
to the zeroth-order terms Pq , Pi = P|, and Pj 12 . 

We will make use of the relations Pj = P£ and 
gx'iX"' _ yx',x"' 

X,X" X:X" 



C X2- 



Otherwise, P* 1 



in equilibrium and P^ = —P$ and 

A X 

x x" = — * '$" f° r ^ e nrs t-o r der correction in V , 
where x IS obtained from \ by the transformation 1 <-> 2. 
For transition from diagonal states in Liouville space to 
diagonal ones we find E*'*, = 0. Finally, we drop all 
E terms which connect states in Liouville space that are 
not compatible, at least to lowest order in Y. As a conse- 
quence the master equations for the linear correction in 
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V read 
= P^;i + P a 1 (Ei;i-E?;i) (Bl) 

= P 2 S 2 ^ 2 + P x I E 1 | 2 — E 2 | 2 J 

+P() ^2 + Pi ( ^l)2 + ^%2 ) + Pi 2 ^12^2 ■ (B2) 

We calculate all diagrams explicitly and find 



yi,i _ fi,i _ _ T 

^1,1 — ^2,2 — Ji 

^2,1 — ^1.1 — ^1,2 ~~ ^2,2 — U 



for V = 0, and 



5,0,1 S i 



1,1 . f>2,l 



J 0,2 



2 



(B3) 
(B4) 



2 ~~ ^2,2 {,12,1 r , .. 

S 12 2 = 7T eV 7 ( e ) Sln 7T ( B5 ) 



2 1Z ' Z 2 

for the first-order correction in V . This leads to 

PI = (B6) 

p2=-^^/'(e)sin|. (B7) 



t]Z']Z = XlZ'lZ =- -iT [1 + /(e)] (B8) 



J ler,ler 



J 2cr,2er 



y^lcr,l(7 y^2(T.lcr ypla.la y^2er,lcr .-, 



J 2o\l<7 



J lcr,2er 



j 2ct,2ct 



(B9) 



for the terms at equilibrium, and 



y,0,lcr _ L 1it,2j + 2j 2er,2cr 



Aio-2o-,i(r r ,. . . 95 

= S l CT 2<x,2<x = 2 eV f ( £ ) Sm 2 

(BIO) 



for the corrections in first order in V. This leads to 



a 1 ; = 



^ = -r v jdj(eW sHv/2y 



(Bll) 

(B12) 



3. Dot-dot interaction, no spin 

Finally we consider the case of a finite intera ctio n be- 
tween dot 1 and 2. For this case, again Eqs. (Bl) and 
(p2|) hold. In the limit /?e ~ 1 but j3U' > I we find 



2. Infinite charging energy 

It is straightforward to generalize the previous discus- 
sion to the case of QDs with infinite charging energy, 
U = oo. In addition to the label for dots 1 and 2, we 
use a to label the spin, and a for the spin opposite to a. 
Since infinite charging energy suppresses all states where 
either one or both dots are occupied with two electrons, 
the HUbert space has nine dimensions. Spin symmetry 
reduces the number of independent diagonal matrix el- 
ements from nine to five, P °, P\° = P^f, P^ = P^ , 
P\%Z = Pilfl, and P^ = Pll%. They are all real. 
Since the transformation 1 <-> 2 and V <-> — V does not 
change the system, the four combinations P®, P^ + P^, 
P\a2a i an< i Pi<t2s are even m ^> while the combination 
P# - P% is odd. 

We are only interested in those off-diagonal matrix el- 
ements for which the projection on the z direction of 
the total spin of the corresponding states is the same. 
There are six of them, but again spin symmetry re- 
duced the number of independent elements to three, 
Pg = PH, and P\°U = P\°l°. Further- 



P1<T 

r 2a 



pl(T 

r 2a i 



yla 



(Pt) 



There are only two independent 
variables, P\% — Pf^ arid Imi^J, which have a linear 
correction in V. 

Again we write P = P + P + ... and £ = £ + £ + 
. . . for the term in zeroth order in V and the first-order 
correction. There are two independent equations which 
look identical to Eqs. (Bl) and (B2) if we replace 1 — > la, 
2 — ► 2cr, and 12 — > lcr2cr, i.e., the equations for different 
spins decouple. 

We calculate all diagrams explicitly and find 



y,l,l y,2,l _ y,l,l 

^2,1 — ^1,1 — ^1 



in equilibrium, and 



-v [i - m] 

«a,i r /3C7' ^ 

2 - E 2 2 = - ln ^— C0S ^ 

7T Z7T 2 



: r 



^0,2 = ^li + ^lil = 2 eV 7'( £ ) sin 2 

^12,2 - u 



(B13) 
(B14) 

(B15) 
(B16) 



for the corrections in first order in V. This leads to 

(B17) 



p nmpuw y p 

1 t l-/(e) 2 2 



^—^^•^sin^/2) (B18) 



with the factor C as defined in Eq. (4.10) 



APPENDIX C: ALTERNATIVE DERIVATION OF 
EQ. (4.7) 

As shown in Ref. p3 the current can be written as 

I=^Jdu,tr{G a r R G r r L }(f L -f R ) (CI) 

In this case for the linear conductance only equilibrium 
Green's functions are involved. These can be determined 
exactly up to all orders, e.g., by using an equations-of- 
motions approach. The result for the retarded Green's 
function is 
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G 



( u-e + iT/2 i(T/2)co S (<p/2)\ ' 
{UJ) ~ \ *(T/2) cos(^/2) u - e + iT/2 J l j 



and the advanced Green's function is the complex conju- 
gate. As a consequence the transmission is 



7» = 



r 2 (^-e) 2 cos 2 (^/2) 

(co-e) 2 + (r/2) 2 (l + cos((^/2)) 2 

1 



(C3) 



(w - e) 2 + (r/2) 2 (l - cos(v?/2)) 2 ' 
An expansion up to first and second order in T yields 

T (1) (uj) = ttTS(uj - e) cos 2 (<^/2) (C4) 



for the first order, which proves Eq. (4.7), and 

T 2 cos 2 (ip/2) 



T (2) (w)=Re 



(w - e + i0+) 2 



(C5) 



for the second order (cotunncling). Note that for ip — 
0, ±27r, ±47r, . . ., the transmission in second order is twice 
the sum of the transmissions through the dots taken 
apart, as expected for constructive interference, while in 
first order fat resonance) no factor 2 is involved. 



1 H. Akera, Phys. Rev. B 47, 6835 (1993). 

2 A. Yacoby, M. Heiblum, D. Mahalu, and H. Shtrikman, 
Phys. Rev. Lett. 74, 4047 (1995). 

3 A.L. Yeyati and M. Biittiker, Phys. Rev. B 52, R14360 
(1995). 

4 G. Hackenbroich and H.A. Weidenmiiller, Phys. Rev. Lett. 
76, 110 (1996). 

5 C. Bruder, R. Fazio, and H. Schoeller, Phys. Rev. Lett. 76, 
114 (1996). 

6 Y. Oreg and Y. Gefen, Phys. Rev. B 55, 13726 (1997). 

7 R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Uman- 
sky, and H. Shtrikman, Nature (London) 385, 417 (1997). 

8 W. Izumida, O. Sakai, and Y. Shimizu, J. Phys. Soc. Jpn. 
66, 717 (1997). 

9 R. Baltin, Y. Gefen, G Hackenbroich, and H.A Wei- 
denmiiller, Eur. Phys. J. B 10, 119 (1999). 

10 R. Baltin and Y. Gefen, Phys. Rev. Lett. 83, 5094 (1999). 

11 Y. Ji, M. Heiblum, D. Sprinzak, D. Mahalu, and H. Shtrik- 
man, Science 290, 779 (2000). 

12 W.G van der Wiel, S. De Franceschi, T. Fujisawa, J.M. 
Elzerman, S. Tarucha, and L.P. Kouwenhoven, Science 
289, 2105 (2000). 

13 D. Loss and E.V. Sukhorukov, Phys. Rev. Lett. 84, 1035 
(2000). 

14 U. Gerland, J. v. Delft, T.A. Costi, and Y. Oreg, Phys. 
Rev. Lett. 84, 3710 (2000). 

15 P.G Silvestrov and Y. Imry, Phys. Rev. Lett. 85, 2565 
(2000). 

16 K. Rang and S.C. Shin, Phys. Rev. Lett. 85, 5619 (2000). 



17 D. Boese, W. Hofstetter, and H. Schoeller, Phys. Rev. B 

64, 125309 (2001). 

18 A.W. Holleitner, C.R. Decker, H. Qin, K. Eberl, and R.H. 
Blick, Phys. Rev. Lett. 87, 25680 (2001). 

19 J. Konig and Y. Gefen, Phys. Rev. Lett. 86, 3855 (2001). 

20 W. Hofstetter, J. Konig, and H. Schoeller, Phys. Rev. Lett. 
87, 156803 (2001). 

21 This is justified so long as the temperature is high enough 
to render the Kondo physics irrelevant. 

22 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 
Abusch-Magder, U. Meirav, M.A. Kastner, Nature (Lon- 
don) 391, 156 (1998). 

23 S.M. Cronenwett, T.H. Oosterkamp, and L.P. Kouwen- 
hoven, Science 281, 540 (1998). 

24 F. Simmel, R.H. Blick, J. P. Kotthaus, W. Wegscheider, and 
M. Bichler, Phys. Rev. Lett. 83, 804 (1999). 

25 J. Schmid, J. Weis, K. Eberl, and K. v. Klitzing, Phys. 
Rev. Lett. 84, 5824 (2000). 

26 L.I. Glazman and M.E. Raikh, Pis'ma Zh. Eksp. Teor. 47, 
378 (1998) [JETP Lett. 47, 452 (1988)]. 

27 T.K. Ng and P.A. Lee, Phys. Rev. Lett. 61, 1768 (1988). 

28 Y. Meir, N.S. Wingreen, and P.A. Lee, Phys. Rev. Lett. 
70, 2601 (1993). 

29 J. Konig, H. Schoeller, and G Schon, Phys. Rev. Lett. 76, 
1715 (1996); J. Konig, J. Schmid, H. Schoeller, and G 
Schon, Phys. Rev. B 54, 16820 (1996). 

30 M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986); L.J. van der 
Pauw, Philips Res. Rep. 13, 1 (1958). 

31 R. Landauer, Philos. Mag. 21, 863 (1970). 

32 M. Biittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. 
Rev. B 31, 6207 (1985); Y. Gefen, Y. Imry, and M.Ya. 
Azbel, Phys. Rev. Lett. 52, 129 (1984). 

33 Y. Meir and N. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 

34 H. Schoeller, in Mesoscopic Electron Transport, edited by 
L.L. Sohn, L.P. Kouwenhoven, and G Schon (Kluwer, 
Dordrecht, 1997); J. Konig, Quantum Fluctuations in the 
Single-Electron Transistor (Shaker, Aachen, 1999). 

35 D.V. Averin and A.A. Odintsov, Phys. Lett. A 140, 251 
(1989); D.V. Averin and Yu.V. Nazarov, Phys. Rev. Lett. 

65, 2446 (1990). 

36 The renormalization terms derived in second order, "co- 
tunneling at resonance" ,lj are negligible in this regime. 

37 J. KSnig, H. Schoeller, and G Schon, Phys. Rev. Lett. 78, 
4482 (1997); Phys. Rev. B 58, 7882 (1998). 

38 B.L. Altshuler, A.G Aronov, and D.E. Khmelnitsky, J. 
Phys. C 15, 7367 (1982). 

39 A. Stern, Y. Aharonov, and Y. Imry, Phys. Rev. A 41, 
3436 (1990). 

40 See also P.A. Mello, Y. Imry, and B. Shapiro, Phys. Rev. 
B 61, 16570 (2000). 

41 The regularization +i0 + is put here by handp-but can be 
derived within a complete cotunneling theoryEZI 

42 Here, we have chosen tiXnte'^ to be positive. If we chose 
a minus sign, all our conclusions would be modified in a 
trivial way, <p — » ip + n. This would modify, though, details 
of the evolution of the transmission phase with the gate 
voltage, which is not considered here. 

43 This can be easily proven within a diagrammatic technique 
as introduced in Appendix H. The diagram of each contri- 
bution to the Green's function G^ ' contains one "double 



21 



vertex" describing direct tunneling through the reference 
arm, say from the left to the right lead. By interchanging 
the labels L and R for the tunnel lines connected to the 
double vertex we can construct another diagram which con- 
tributes to G^ 1 ' as well. The double vertex now describes 
tunneling from right to left. Since in equilibrium the left 
and right leads are described by the same Fermi function, 
the analytic expressions for the new diagram are almost 
identical to the old ones. The only difference is a reversal 
of the AB phase, ip — > —ip, due to the exchange of the dou- 
ble vertex. This proves the desired symmetry relation. It 
is straightforward to generalize this proof for Green's func- 
tion G^ n ' of arbitrary high order in t Te and, thus, for the 
total Green's function as well. 

B.L. Altshuler, Y. Gefen, A. Kamenev, and L.S. Levitov, 
Phys. Rev. Lett. 78, 2803 (1997). 

T.V. Shahbazyan and M.E. Raikh, Phys. Rev. B 49, 17123 
(1994). 

' J. Konig, Y. Gefen, and G. Schon, Phys. Rev. Lett. 81, 
4468 (1998). 

We consider here the generic case where the signs of all tun- 
neling matrix elements ti,,dati, th,dot2, tfl,dotl, and tn,dot2 
are the same (in the absence of flux). Another choice of 
the relative phases among them will modify the detailed 
conclusions of the following analysis. 

We thank Yang Ji and M. Heiblum for discussions on this 
point. 

' For completeness we mention that further minus signs can 
appear in diagrams which include states with two or more 
fermions. They are due to the fact that a two-vertex con- 
tribution such as CfCJ| |) leads to a minus sign, — | J.) (see 
Refs.Eslfjl). 



22 



